VapC10 toxin of the legume symbiont Sinorhizobium meliloti targets tRNASer and controls intracellular lifestyle.

The soil bacterium Sinorhizobium meliloti can establish a nitrogen-fixing symbiosis with the model legume Medicago truncatula. The rhizobia induce the formation of a specialized root organ called nodule, where they differentiate into bacteroids and reduce atmospheric nitrogen into ammonia. Little is known on the mechanisms involved in nodule senescence onset and in bacteroid survival inside the infected plant cells. Although toxin-antitoxin (TA) systems have been shown to promote intracellular survival within host cells in human pathogenic bacteria, their role in symbiotic bacteria was rarely investigated. S. meliloti encodes several TA systems, mainly of the VapBC family. Here we present the functional characterization, through a multidisciplinary approach, of the VapBC10 TA system of S. meliloti. Following a mapping by overexpression of an RNase in Escherichia coli (MORE) RNA-seq analysis, we demonstrated that the VapC10 toxin is an RNase that cleaves the anticodon loop of two tRNASer. Thereafter, a bioinformatics approach was used to predict VapC10 targets in bacteroids. This analysis suggests that toxin activation triggers a specific proteome reprogramming that could limit nitrogen fixation capability and viability of bacteroids. Accordingly, a vapC10 mutant induces a delayed senescence in nodules, associated to an enhanced bacteroid survival. VapBC10 TA system could contribute to S. meliloti adaptation to symbiotic lifestyle, in response to plant nitrogen status.

Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

AvP: A software package for automatic phylogenetic detection of candidate horizontal gene transfers.

Horizontal gene transfer (HGT) is the transfer of genes between species outside the transmission from parent to offspring. Due to their impact on the genome and biology of various species, HGTs have gained broader attention, but high-throughput methods to robustly identify them are lacking. One rapid method to identify HGT candidates is to calculate the difference in similarity between the most similar gene in closely related species and the most similar gene in distantly related species. Although metrics on similarity associated with taxonomic information can rapidly detect putative HGTs, these methods are hampered by false positives that are difficult to track. Furthermore, they do not inform on the evolutionary trajectory and events such as duplications. Hence, phylogenetic analysis is necessary to confirm HGT candidates and provide a more comprehensive view of their origin and evolutionary history. However, phylogenetic reconstruction requires several time-consuming manual steps to retrieve the homologous sequences, produce a multiple alignment, construct the phylogeny and analyze the topology to assess whether it supports the HGT hypothesis. Here, we present AvP which automatically performs all these steps and detects candidate HGTs within a phylogenetic framework.

Improvement of eukaryotic protein predictions from soil metagenomes.

During the last decades, metagenomics has highlighted the diversity of microorganisms from environmental or host-associated samples. Most metagenomics public repositories use annotation pipelines tailored for prokaryotes regardless of the taxonomic origin of contigs. Consequently, eukaryotic contigs with intrinsically different gene features, are not optimally annotated. Using a bioinformatics pipeline, we have filtered 7.9 billion contigs from 6,872 soil metagenomes in the JGI's IMG/M database to identify eukaryotic contigs. We have re-annotated genes using eukaryote-tailored methods, yielding 8 million eukaryotic proteins and over 300,000 orphan proteins lacking homology in public databases. Comparing the gene predictions we made with initial JGI ones on the same contigs, we confirmed our pipeline improves eukaryotic proteins completeness and contiguity in soil metagenomes. The improved quality of eukaryotic proteins combined with a more comprehensive assignment method yielded more reliable taxonomic annotation. This dataset of eukaryotic soil proteins with improved completeness, quality and taxonomic annotation reliability is of interest for any scientist aiming at studying the composition, biological functions and gene flux in soil communities involving eukaryotes.

Movements of transposable elements contribute to the genomic plasticity and species diversification in an asexually reproducing nematode pest.

Despite reproducing without sexual recombination, Meloidogyne incognita is an adaptive and versatile phytoparasitic nematode. This species displays a global distribution, can parasitize a large range of plants, and can overcome plant resistance in a few generations. The mechanisms underlying this adaptability remain poorly known. At the whole-genome level, only a few single nucleotide variations have been observed across different geographical isolates with distinct ranges of compatible hosts. Exploring other factors possibly involved in genomic plasticity is thus important. Transposable elements (TEs), by their repetitive nature and mobility, can passively and actively impact the genome dynamics. This is particularly expected in polyploid hybrid genomes such as the one of M. incognita. Here, we have annotated the TE content of M. incognita, analyzed the statistical properties of this TE landscape, and used whole-genome pool-seq data to estimate the mobility of these TEs across twelve geographical isolates, presenting variations in ranges of compatible host plants. DNA transposons are more abundant than retrotransposons, and the high similarity of TE copies to their consensus sequences suggests they have been at least recently active. We have identified loci in the genome where the frequencies of presence of a TE showed substantial variations across the different isolates. Overall, variations in TE frequencies across isolates followed their phylogenetic divergence, suggesting TEs participate in the species diversification. Compared with the M. incognita reference genome, we detected isolate and lineage-specific de novo insertion of some TEs, including within genic regions or in the upstream regulatory regions. We validated by PCR the insertion of some of these TEs inside genic regions, confirming TE movements have possible functional impacts. Overall, we show DNA transposons can drive genomic plasticity in M. incognita and their role in genome evolution of other parthenogenetic animal deserves further investigation.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes.

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.

Sequencing the Genome of Indian Flying Fox, Natural Reservoir of Nipah Virus, Using Hybrid Assembly and Conservative Secondary Scaffolding.

Indian fruit bats, flying fox Pteropus medius was identified as an asymptomatic natural host of recently emerged Nipah virus, which is known to induce a severe infectious disease in humans. The absence of P. medius genome sequence presents an important obstacle for further studies of virus-host interactions and better understanding of mechanisms of zoonotic viral emergence. Generation of the high-quality genome sequence is often linked to a considerable effort associated to elevated costs. Although secondary scaffolding methods have reduced sequencing expenses, they imply the development of new tools for the integration of different data sources to achieve more reliable sequencing results. We initially sequenced the P. medius genome using the combination of Illumina paired-end and Nanopore sequencing, with a depth of 57.4x and 6.1x, respectively. Then, we introduced the novel scaff2link software to integrate multiple sources of information for secondary scaffolding, allowing to remove the association with discordant information among two sources. Different quality metrics were next produced to validate the benefits from secondary scaffolding. The P. medius genome, assembled by this method, has a length of 1,985 Mb and consists of 33,613 contigs and 16,113 scaffolds with an NG50 of 19 Mb. At least 22.5% of the assembled sequences is covered by interspersed repeats already described in other species and 19,823 coding genes are annotated. Phylogenetic analysis demonstrated the clustering of P. medius genome with two other Pteropus bat species, P. alecto and P. vampyrus, for which genome sequences are currently available. SARS-CoV entry receptor ACE2 sequence of P. medius was 82.7% identical with ACE2 of Rhinolophus sinicus bats, thought to be the natural host of SARS-CoV. Altogether, our results confirm that a lower depth of sequencing is enough to obtain a valuable genome sequence, using secondary scaffolding approaches and demonstrate the benefits of the scaff2link application. The genome sequence is now available to the scientific community to (i) proceed with further genomic analysis of P. medius, (ii) to characterize the underlying mechanism allowing Nipah virus maintenance and perpetuation in its bat host, and (iii) to monitor their evolutionary pathways toward a better understanding of bats' ability to control viral infections.

Social learning in great white pelicans (Pelecanus onocrotalus): A preliminary study.

Great white pelicans (Pelecanus onocrotalus) exhibit life-history parameters and ecological traits thought to be associated with social learning, and advanced cognitive processing more generally. In this study we investigated whether this species can acquire novel behavior socially in a foraging context. Birds from the test group watched a trained conspecific opening an opaque box containing a food reward by using its beak, whereas the control group had no demonstrator but saw the box for an equivalent time span. Individuals from both groups were subsequently allowed access to the box. Subjects of the test group performed significantly better than the control group. This is the first experimental evidence of social learning in a cooperatively hunting bird. Further studies are needed in order to shed light on the factors favoring the evolution of this capacity, by testing different pelican species that vary in their ecology.

Population genomics supports clonal reproduction and multiple independent gains and losses of parasitic abilities in the most devastating nematode pest.

The root-knot nematodes are the most devastating worms to worldwide agriculture with Meloidogyne incognita being the most widely distributed and damaging species. This parasitic and ecological success seems surprising given its supposed obligatory clonal reproduction. Clonal reproduction has been suspected based on cytological observations but, so far, never confirmed by population genomics data. As a species, M. incognita is highly polyphagous with thousands of host plants. However, different M. incognita isolates present distinct and overlapping patterns of host compatibilities. Historically, four "host races" had been defined as a function of ranges of compatible and incompatible plants. In this study, we used population genomics to assess whether (a) reproduction is actually clonal in this species, (b) the host races follow an underlying phylogenetic signal or, rather represent multiple independent transitions, and (c) how genome variations associate with other important biological traits such as the affected crops and geographical distribution. We sequenced the genomes of 11 M. incognita isolates across Brazil that covered the four host races in replicates. By aligning the genomic reads of these isolates to the M. incognita reference genome assembly, we identified point variations. Analysis of linkage disequilibrium and 4-gametes test showed no evidence for recombination, corroborating the clonal reproduction of M. incognita. The few point variations between the isolates showed no significant association with the host races, the geographical origin of the samples, or the crop on which they have been collected. Addition of isolates from other locations around the world confirmed this lack of underlying phylogenetic signal. This suggests multiple gains and losses of parasitic abilities and adaptations to different environments account for the broad host spectrum and wide geographical distribution of M. incognita and thus to its high economic impact. This surprising adaptability without sex poses both evolutionary and agro-economic challenges.

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic, plant-parasitic nematode Meloidogyne incognita.

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root-knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant-parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.

The Landscape of L1 Retrotransposons in the Human Genome Is Shaped by Pre-insertion Sequence Biases and Post-insertion Selection.

L1 retrotransposons are transposable elements and major contributors of genetic variation in humans. Where L1 integrates into the genome can directly impact human evolution and disease. Here, we experimentally induced L1 retrotransposition in cells and mapped integration sites at nucleotide resolution. At local scales, L1 integration is mostly restricted by genome sequence biases and the specificity of the L1 machinery. At regional scales, L1 shows a broad capacity for integration into all chromatin states, in contrast to other known mobile genetic elements. However, integration is influenced by the replication timing of target regions, suggesting a link to host DNA replication. The distribution of new L1 integrations differs from those of preexisting L1 copies, which are significantly reshaped by natural selection. Our findings reveal that the L1 machinery has evolved to efficiently target all genomic regions and underline a predominant role for post-integrative processes on the distribution of endogenous L1 elements.

Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes.

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

How Long Does Wolbachia Remain on Board?

Wolbachia bacteria infect about half of all arthropods, with diverse and extreme consequences ranging from sex-ratio distortion and mating incompatibilities to protection against viruses. These phenotypic effects, combined with efficient vertical transmission from mothers to offspring, satisfactorily explain the invasion dynamics of Wolbachia within species. However, beyond the species level, the lack of congruence between the host and symbiont phylogenetic trees indicates that Wolbachia horizontal transfers and extinctions do happen and underlie its global distribution. But how often do they occur? And has the Wolbachia pandemic reached its equilibrium? Here, we address these questions by inferring recent acquisition/loss events from the distribution of Wolbachia lineages across the mitochondrial DNA tree of 3,600 arthropod specimens, spanning 1,100 species from Tahiti and surrounding islands. We show that most events occurred within the last million years, but are likely attributable to individual level variation (e.g., imperfect maternal transmission) rather than population level variation (e.g., Wolbachia extinction). At the population level, we estimate that mitochondria typically accumulate 4.7% substitutions per site during an infected episode, and 7.1% substitutions per site during the uninfected phase. Using a Bayesian time calibration of the mitochondrial tree, these numbers translate into infected and uninfected phases of approximately 7 and 9 million years. Infected species thus lose Wolbachia slightly more often than uninfected species acquire it, supporting the view that its present incidence, estimated here slightly below 0.5, represents an epidemiological equilibrium.

Full-field interferometry for counting and differentiating aquatic biotic nanoparticles: from laboratory to Tara Oceans.

There is a huge abundance of viruses and membrane vesicles in seawater. We describe a new full-field, incoherently illuminated, shot-noise limited, common-path interferometric detection method that we couple with the analysis of Brownian motion to detect, quantify, and differentiate biotic nanoparticles. We validated the method with calibrated nanoparticles and homogeneous DNA or RNA viruses. The smallest virus size that we characterized with a suitable signal-to-noise ratio was around 30 nm in diameter. Analysis of Brownian motions revealed anisotropic trajectories for myoviruses.We further applied the method for vesicles detection and for analysis of coastal and oligotrophic samples from Tara Oceans circumnavigation.

SENCA: A Multilayered Codon Model to Study the Origins and Dynamics of Codon Usage.

Gene sequences are the target of evolution operating at different levels, including the nucleotide, codon, and amino acid levels. Disentangling the impact of those different levels on gene sequences requires developing a probabilistic model with three layers. Here we present SENCA (site evolution of nucleotides, codons, and amino acids), a codon substitution model that separately describes 1) nucleotide processes which apply on all sites of a sequence such as the mutational bias, 2) preferences between synonymous codons, and 3) preferences among amino acids. We argue that most synonymous substitutions are not neutral and that SENCA provides more accurate estimates of selection compared with more classical codon sequence models. We study the forces that drive the genomic content evolution, intraspecifically in the core genome of 21 prokaryotes and interspecifically for five Enterobacteria. We retrieve the existence of a universal mutational bias toward AT, and that taking into account selection on synonymous codon usage has consequences on the measurement of selection on nonsynonymous substitutions. We also confirm that codon usage bias is mostly driven by selection on preferred codons. We propose new summary statistics to measure the relative importance of the different evolutionary processes acting on sequences.

Large-scale nucleotide optimization of simian immunodeficiency virus reduces its capacity to stimulate type I interferon in vitro.

UNLABELLED: Lentiviral RNA genomes present a strong bias in their nucleotide composition with extremely high frequencies of A nucleotide in human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV). Based on the observation that human optimization of RNA virus gene fragments may abolish their ability to stimulate the type I interferon (IFN-I) response, we identified the most biased sequences along the SIV genome and showed that they are the most potent IFN-I stimulators. With the aim of designing an attenuated SIV genome based on a reduced capacity to activate the IFN-I response, we synthesized artificial SIV genomes whose biased sequences were optimized toward macaque average nucleotide composition without altering their regulatory elements or amino acid sequences. A synthetic SIV optimized with 169 synonymous mutations in gag and pol genes showed a 100-fold decrease in replicative capacity. Interestingly, a synthetic SIV optimized with 70 synonymous mutations in pol had a normal replicative capacity. Its ability to stimulate IFN-I was reduced when infected cells were cocultured with reporter cells. IFN regulatory factor 3 (IRF3) transcription factor was required for IFN-I stimulation, implicating cytosolic sensors in the detection of SIV-biased RNA in infected cells. No reversion of introduced mutations was observed for either of the optimized viruses after 10 serial passages. In conclusion, we have designed large-scale nucleotide-modified SIVs that may display attenuated pathogenic potential. IMPORTANCE: In this study, we synthesized artificial SIV genomes in which the most hyperbiased sequences were optimized to bring them closer to the nucleotide composition of the macaque SIV host. Interestingly, we generated a stable synthetic SIV optimized with 70 synonymous mutations in pol gene, which had a normal replicative capacity but a reduced ability to stimulate type I IFN. This demonstrates the possibility to rationally change viral nucleotide composition to design replicative and genetically stable lentiviruses with attenuated pathogenic potentials.

The biased nucleotide composition of HIV-1 triggers type I interferon response and correlates with subtype D increased pathogenicity.

The genome of human immunodeficiency virus (HIV) has an average nucleotide composition strongly biased as compared to the human genome. The consequence of such nucleotide composition on HIV pathogenicity has not been investigated yet. To address this question, we analyzed the role of nucleotide bias of HIV-derived nucleic acids in stimulating type-I interferon response in vitro. We found that the biased nucleotide composition of HIV is detected in human cells as compared to humanized sequences, and triggers a strong innate immune response, suggesting the existence of cellular immune mechanisms able to discriminate RNA sequences according to their nucleotide composition or to detect specific secondary structures or linear motifs within biased RNA sequences. We then extended our analysis to the entire genome scale by testing more than 1300 HIV-1 complete genomes to look for an association between nucleotide composition of HIV-1 group M subtypes and their pathogenicity. We found that subtype D, which has an increased pathogenicity compared to the other subtypes, has the most divergent nucleotide composition relative to the human genome. These data support the hypothesis that the biased nucleotide composition of HIV-1 may be related to its pathogenicity.

Mammalian X chromosome inactivation evolved as a dosage-compensation mechanism for dosage-sensitive genes on the X chromosome.

How and why female somatic X-chromosome inactivation (XCI) evolved in mammals remains poorly understood. It has been proposed that XCI is a dosage-compensation mechanism that evolved to equalize expression levels of X-linked genes in females (2X) and males (1X), with a prior twofold increase in expression of X-linked genes in both sexes ("Ohno's hypothesis"). Whereas the parity of X chromosome expression between the sexes has been clearly demonstrated, tests for the doubling of expression levels globally along the X chromosome have returned contradictory results. However, changes in gene dosage during sex-chromosome evolution are not expected to impact on all genes equally, and should have greater consequences for dosage-sensitive genes. We show that, for genes encoding components of large protein complexes (≥ 7 members)--a class of genes that is expected to be dosage-sensitive--expression of X-linked genes is similar to that of autosomal genes within the complex. These data support Ohno's hypothesis that XCI acts as a dosage-compensation mechanism, and allow us to refine Ohno's model of XCI evolution. We also explore the contribution of dosage-sensitive genes to X aneuploidy phenotypes in humans, such as Turner (X0) and Klinefelter (XXY) syndromes. X aneuploidy in humans is common and is known to have mild effects because most of the supernumerary X genes are inactivated and not affected by aneuploidy. Only genes escaping XCI experience dosage changes in X-aneuploidy patients. We combined data on dosage sensitivity and XCI to compute a list of candidate genes for X-aneuploidy syndromes.

An externally modulated, noise-driven switch for the regulation of SPI1 in Salmonella enterica serovar Typhimurium.

In this work we consider the regulation system present on the SPI1 pathogenicity island of Salmonella enterica serovar Typhimurium. It is well-known that HilA is the central regulator in the overall scheme of SPI1 regulation and directly binds to virulence operons and activates their expression. The regulation of the expression of HilA is via a complex feed-forward loop involving three transcriptional activators: HilC, HilD and RtsA, and the negative regulator HilE. Our aim is to model this regulation network and study its dynamical behavior. We show that this regulatory system can display a bistable behavior relevant to the biology of Salmonella, and that noise can be a driving force in this system.

Inference of sparse combinatorial-control networks from gene-expression data: a message passing approach.

BACKGROUND: Transcriptional gene regulation is one of the most important mechanisms in controlling many essential cellular processes, including cell development, cell-cycle control, and the cellular response to variations in environmental conditions. Genes are regulated by transcription factors and other genes/proteins via a complex interconnection network. Such regulatory links may be predicted using microarray expression data, but most regulation models suppose transcription factor independence, which leads to spurious links when many genes have highly correlated expression levels. RESULTS: We propose a new algorithm to infer combinatorial control networks from gene-expression data. Based on a simple model of combinatorial gene regulation, it includes a message-passing approach which avoids explicit sampling over putative gene-regulatory networks. This algorithm is shown to recover the structure of a simple artificial cell-cycle network model for baker's yeast. It is then applied to a large-scale yeast gene expression dataset in order to identify combinatorial regulations, and to a data set of direct medical interest, namely the Pleiotropic Drug Resistance (PDR) network. CONCLUSIONS: The algorithm we designed is able to recover biologically meaningful interactions, as shown by recent experimental results 1. Moreover, new cases of combinatorial control are predicted, showing how simple models taking this phenomenon into account can lead to informative predictions and allow to extract more putative regulatory interactions from microarray databases.